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Abstract 

We give an overview of some properties of Lie algebras generated by at most 5 extremal 
elements. In particular, for any finite graph T and any field K of characteristic not 2, we 
consider an algebraic variety X over K whose if-points parametrize Lie algebras generated 
by extremal elements. Here the generators correspond to the vertices of the graph, and we 
prescribe commutation relations corresponding to the nonedges of T. 

We show that, for all connected undirected finite graphs on at most 5 vertices, X is a finite- 
dimensional affine space. Furthermore, we show that for maximal-dimensional Lie algebras 
generated by 5 extremal elements, X is a single point. The latter result implies that the 
bilinear map describing extremality must be identically zero, so that all extremal elements 
are sandwich elements and the only Lie algebra of this dimension that occurs is nilpotent. 

These results were obtained by extensive computations with the Magma computational 
algebra system. The algorithms developed can be applied to arbitrary T (i.e., without restric- 
tion on the number of vertices), and may be of independent interest. 

1 Introduction 

We assume throughout the paper that K is a field of characteristic distinct from 2. Let L be a 
Lie algebra over K. A non-zero element x G L is called extremal if [x, [x, y]] € Kx for all y € L. if 
x is extremal, the existence of a linear map f x :L^K such that [x, [x, y]] = f x {y)x for all y e L 
immediately follows from linearity of [■,■]. If for some extremal element x € L this linear map f x 
is identically 0, we call x a sandwich. 

Extremal elements were originally introduced by Chernousov [Che89] in his proof of the Hassc 
principle for Eg. Zcl'manov and Kostrikin proved that, for every n, the universal Lie algebra L n 
generated by a finite number of sandwich elements X\, . . . , x n is finite-dimensional [ZK90]. Cohen, 
Steinbach, Ushirobira, and Wales generalized this result and proved that a Lie algebra generated 
by a finite number of extremal elements is finite dimensional. Moreover, they give an explicit lower 
bound on the number of extremal elements required to generate each of the classical Lie algebras 
[CSUW01]. Recently, in 't panhuis, Postma, and the author gave explicit presentations for Lie 
algebras of type A„, B„, C„, and D„, by means of minimal sets of extremal generators [itpPR09]. 

Moreover, Draisma and in 't panhuis considered finite graphs and corresponding algebraic 
varieties whose points parametrize Lie algebras generated by extremal elements. They proved in 
particular that if the graph is a simply laced Dynkin diagram of affine type, all points in an open 
dense subset of the affine variety parametrize Lie algebras isomorphic to the simple Chevalley Lie 
algebra corresponding to the associated Dynkin diagram of finite type [Ditp08] . 

Looking at these Lie algebras from a different point of view, Cohen, Ivanyos, and the author 
proved that if L is a Lie algebra over a field K (of characteristic distinct from 2 and 3) that has 
an extremal element that is not a sandwich, then L is generated by extremal elements, with one 
exception in characteristic 5 [CIR08]. 



The strong connection between extremal elements and geometries is further investigated in two 
papers by Cohen and Ivanyos [CI06, CI07], and in the Ph.D. theses by Postma [Pos07] and in 't 
panhuis [itp09]. 

2 Preliminaries and main results 

We follow the setup of [itpPR09] and [Ditp08]. Assume that T is a connected undirected finite 
graph with n vertices, without loops or multiple bonds, and that if is a field of characteristic 
distinct from 2. We let II be the vertex set of T and denote adjacency of two vertices x, y £ II by 
i ~ j. 

We denote by T(K, T) (often abbreviated to T) the quotient of the free Lie algebra over K 
generated by II modulo the relations 

[x, y] = for all x,y £U with i/y. 

Often, we write elements of T as linear combinations of monomials [x\, [x 2 , . . . , [x/_i,x;]]], where 
xi, . . . , xi £ II. The length of such a monomial is said to be I, and we often abbreviate such a 
monomial to [xi, . . . ,x{\. Note that T inherits the natural N-grading from the free Lie algebra 
generated by II; homogeneous elements of T are linear combinations of monomials of equal length. 

We write F* for the space of all i-T-linear functions F — > K. For every / £ (J 7 *) 11 , also written 
(fx)xen, w e denote by C(K, T, /) (often abbreviated to £(/)) the quotient of T(K, V) by the ideal 
1(f) generated by the infinitely many elements 

[x, [x, y}} - f x (y)x for x e II and y £ T. 

By construction C(f ) is a Lie algebra generated by |IT| = n extremal elements, the extremal 
generators corresponding to the vertices of T and commuting whenever they are not adjacent. The 
element f x £ J-* is a parameter expressing the extremality of x £ II. 

In the Lie algebra C(0) the elements of II map to sandwich elements. By [ZK90] this Lie algebra 
is finite-dimensional; for general / £ (J 7 *) 11 we have dim(£(/)) < dim(£(0)) by [CSUW01, Lemma 
4.3]. It is therefore natural to focus on the Lie algebras C(f) of maximal possible dimension, 
i.e., those of dimension dim(£(0)). We define the set 

X := {/ e (F*) U I dim(£(/)) = dim(£(0))}, 

the parameter space for all maximal-dimensional Lie algebras of the form £(/). 

Example 1. Consider the case where T consists of two vertices x, y, connected by an edge. 
Then for every / = (f x , f y ) £ (F*) n the Lie algebra C(f) is spanned by B — {x, y, [x, y]} since 
[x, [x,y]] = f x (y)x (where f x (y) £ K), and similarly [y, [x,y]] = -f y (x)y. 

In particular, for / = (0,0) we find that C(f) is the Heisenberg algebra, which is 3- 
dimensional. For general /, the requirement that C(f) is 3-dimensional implies (using the 
Jacobi identity) 

f y (x)[x,y] = [x, [y, [y,x]]] = [y, [x, [y,x]}] + [[x,y], [y,x\] = 

- [y, [x, [x,y]]} + = -f x (y)[y,x] = f x (y)[y,x], 

so that f x (y) must be equal to f y (x) since we assumed [x,y] =/= 0. 

Consequently, 3-dimensional Lie algebras generated by a distinguished pair of extremal gen- 
erators are parametrized by the single value f x (y). As mentioned before, if f x (y) = then 
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£(/) is the Heisenberg algebra. It is straightforward to verify that those Lie algebras where 
f x (y) 7^ are mutually isomorphic and isomorphic to the split simple Lie algebra of type Ai. 
The parameter space A defined above, then, is the affine line, and all Lie algebras corresponding 
to the nonzero points on the line are mutually isomorphic. 

The following theorem asserts that the two generator case is exemplary: A always carries the 
structure of an affine algebraic variety. 

Theorem 2 ([Ditp08, Theorem 1]). The set X is naturally the set of K -rational points of an 
affine variety of finite type defined over K . This variety can be described as follows. Fix any 
finite- dimensional subspace V of T consisting of homogeneous elements such that V +1(0) = T . 
Then the restriction map 

X (V*) U , f ^ (/* |v)xen 
maps X injectively onto the set of K -rational points of a closed subvariety of (V*) n . This yields 
a K-variety structure of X which is independent of the choice of V . 

Example 1 (continued). We give an example of an embedding of X as mentioned in 
Theorem 2. Recall that we showed X to be the affine line; we may take V = B = {x, y, [x, y]}. 
Then, for / £ X we have 

/ !-> ((fx(x), f x (y), f x ([x,y})), (f y (x), f y (y), f y ([x,y])), 

and if we again identify elements / £ X with f x (y) € K we find / >-> ((0, /, 0), (/, 0, 0)). 

We have now introduced enough background to be able to state our main results. 

Theorem 3. Let T be the complete graph on 5 vertices, K a field of characteristic not 2, and, 
as above, the set X the parameter space for all maximal-dimensional Lie algebras of the form 
C(K,T,f). Then X = {0}. 

The following corollary immediately follows from this theorem and [CSUW01, Lemma 4.2] 
(and Section 6 with regards to the dimensions mentioned). 

Corollary 4. Suppose L is a Lie algebra generated by 5 extremal elements over the field K , where 
char(A) 7^ 2. If L is of maximal dimension among such Lie algebras (i.e., dim(i) = 537 if 
char(A) ^ 3 and dim(L) = 538 i/char(A) = 3) then L is nilpotent. 

Looking at results for other cases shows why this result is interesting. For example, if F is 
the complete graph on on 2, 3, or 4 vertices, then X is nontrivial. Moreover, for almost all 
/ £ X, £(/) is a simple Lie algebra of type A l7 A 2 , and D 4 , respectively (see [CSUW01, Roo05] 
and Section 6 of this paper). Moreover, there are infinite families of graphs T, for example the 
affine Dynkin diagrams, where C(f) is almost always a simple Lie algebra of classical type (see 
[itpPR09, Ditp08]). In short: this is the first nontrivial case encountered where X is a point. 

Note that the fact that A is a point means that the extremality of the generators and the 
maximal-dimensionality assumption force all generators to be sandwiches, and therefore force L 
to be nilpotent. 

We mention one more result. 

Theorem 5. Let T be a connected undirected finite graph with at most 5 vertices, without loops 
or multiple bonds, K a field of characteristic not 2, and, as above, the set X the parameter space 
for all maximal- dimensional Lie algebras of the form C(K,T, f ). Then X is isomorphic to a 
finite- dimensional affine space over K. The dimension of X is given in Tables 1-3. 
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This theorem provides a partial answer to one of the questions posed in [Ditp08, Section 5.2], 
namely, "Is X always an affine space?" . Similarly, the results discussed in Section 6 suggest that in 
the cases we considered one of the other questions posed, "Is there always a generic Lie algebra?" 
can be answered affirmatively. 

The remainder of this paper 

In Section 3 we present a number of known results on these Lie algebras and their elements that 
we need for our proofs and calculations. In Section 4 we introduce /-sets, a concept we need in 
order to be able to calculate with the affine algebraic variety X that parametrizes the Lie algebras 
we are interested in. In Section 5 wc discuss the algorithms we developed and implemented in 
Magma to obtain the computation results presented in Section 6. Finally, in Section 7 we briefly 
reflect on our results and possible future research. 

3 Some properties of extremal elements 

In this section we present a number of properties of extremal elements and the Lie algebras 
generated by them. The following identities go back to Premet and are commonly called the 
Premet identities. We have, for x an extremal element and for all j/,z£l: 



There is a number of properties of the functions f x introduced earlier that we should discuss. 
The following result is an important tool in our computations. 

Theorem 8 ([CSUW01, Theorem 2.5]). Suppose that L is a Lie algebra over K generated by 
extremal elements. There is a unique bilinear symmetric form f : L x L — > K such that, for each 
extremal element x € L, the linear form f x coincides with y f(x,y). This form is associative, 
in the sense that f(x, [y, z}) = f([x, y], z) for all x,y,z e L. 

This theorem implies a number of identities in the values of /, as presented in the following 
lemma. We use these relations extensively in our computer calculations. 

Lemma 9. Let x, y, z,y\, . . . ,yi <G II, let (f x ) xe n € (J 7 *) 11 , and let m e £(/). Then the following 
equalities hold: 



(i) fx(y) — whenever x — y or x ^ y, 

(H) fx{y) = fy(x), 

(Hi) f x ([y,m]) — whenever x — y or x ^ y, 

(iv) fx{[y,m]) = -f y (x,m), 

(v) f x {[y, [z,x]}) = f x {y)f x (z), 

(vi) fx([y, [x,m]]) = -f x {y)f x {m), 

(vii) f x ([y, [z,m]]) = f z ([y, [x,m]]) - f z ([x, [y,m]}), 

(via) f x {[yi, [y 2 , [yi-i,yi]}]) = (-l)' -1 /^ ([yi-u [yi-2, [yi, [yi,x]]}]). 



2[x, [y, [x,z]}} = f x ([y,z])x - f x (z)[x,y] - f x (y)[x,z], 
2[[x,y], [x,z]j = f x ([y,z])x + f x (z)[x,y] - f x (y)[x,z}. 



(6) 
(7) 
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Proof (i) follows from the observation that [x, [x,y]] = whenever x = y or x ^ y. (ii) follows 
immediately from the symmetry of / in Theorem 8: f x {y) = f{x,y) = f{y,x) = f y (x). (hi) 
and (iv) follow similarly from its associativity: f x ([y,m]) = f(x,[y 7 m]) — f([x,y],m), which 
is equal to for case (hi) and reduces to —f([y, x], to) = — f(y, [x, to]) = — f y {[x,m]) in the 
case of (iv). For (v), observe f x ([y, [z,x]}) = f(x, [y, [z,x]}) = -f(y, [x, [z,x]}) = f(y,f x (z)x) = 
fx{z)f{y 1 x) — fx(z),f x {y) by associativity and bilinearity of /; (vi) follows similarly, (vii) follows 
from the Jacobi identity and associativity: f x {[y 1 [z, to]]) = f(x, [y, [z,m]]) — — f([x,y], [to, z]) = 
-f([[x,y],m],z) = -f z ([[x,y],m]) = -f z ([x, [y,m]]) + f z ([y, [a;, to]]). Finally, to see (viii) simply 
apply associativity of / and anti-symmetry of Lie algebras I — 1 times. □ 

4 /-sets 

Our goal in the research described here was to find more information on the structure of the 
variety X describing the parameter space, and of £(/) for various / el. In order to achieve this 
goal, we have to compute a multiplication table for £(/). There are two ways to approach this 
problem: either pick an / e X in advance and compute a multiplication table for £(/), or compute 
a general multiplication table, i.e., one that has entries in the coordinate ring of X. The latter 
approach, which we have chosen to pursue, has two important advantages: Firstly, we only need 
to compute the multiplication table once, and can afterwards easily instantiate it for particular 
/ € X. Secondly, by insisting that C(f) be a Lie algebra and of maximal dimension we obtain 
information about the structure of X. 

Example 1 (continued). In this example we would automatically recover the relation 
fx(y) = fy(x) via the Jacobi identity on y, x, and [x, y]: 

= [y, [x, [x,y]}} + [x, [[x,y],y}\ + [[x,y], [y,x]} = f x (y)[y,x] + f y (x)[x,y] + 0, 

and f x ([x,y]) — via evaluating [ y]]] from the inside out: 

f x ([x,y])x = [x, [x, [x,y]}} = [x,f x (y)x] = f x (y)[x,x] = 0, 

so that f x ([x, y]) = since the fact that £(/) is of maximal dimension implies that x ^ 0. 

To aid in the description of our algorithms we introduce the concept of f-sets. For the re- 
mainder of this section fix a graph T and a field K, and let B be a set of monomials of T(K, T) 
that project to a monomial if-basis of £(0) (which exists by [itpPR09, Lemma 3.1]). Further- 
more, let F be some subset of II x B, and let Rp be the rank |F| multivariate polynomial ring 
K[f y (c) | (y,c) e F] with variables f y (c). Furthermore, let r : LT x B — > Rp, (x,b) r x (b) be a 
map from LI x B into this polynomial ring. We call such a pair (F,r) an f-set of size |F|. 

Moreover, if r is such that for all / £ X we have, under the evaluation f y (c) f y (c) for 

(y,c)eF, 

r x (b) = f x (b) for all x e II, b e B, 

then we call (F, r) a sufficient f-set. 

Finally, suppose (F,r) is an /-set where we require, to exclude trivialities, that r y (c) = f y (c) 
for all (y, c) e F. We construct, for an arbitrary v £ K F , an element (f x ) xe n of (J 7 *) 11 as follows. 
Take, for every b <G B, the value f x (b) £ K to be r x (b) € Rf evaluated in v, that is, f y (c) i— > V( y c y 
This defines a map ip : K F — > (J 7 *) 11 . Note that £(ip(v)) is always a Lie algebra, but not always 
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one of maximal dimension. If, however, £(ip(v)) is a Lie algebra of dimension dim(£(0)) for all 
v e K F , then we call (F, r) a free f-set. 

For example, (II x B, (y, c) i-> fy(c)) is trivially a sufficient f-set. It is, however, in general not 
a free /-set. On the other hand, (0, (y, c) ^ 0) forms a free /-set, but in general not a sufficient 
one. 

Example 1 (continued). Recall LT = {x,y} and B = {x,y,[x,y]}, so that II x B = 
{(x, x), (x, y), {x, [x, y]), (y, x), (y, y), (y, [x, y})}. Then (II x B, (y, c) ^ t y (c)) is a (trivial) suffi- 
cient /-set of size 6. In order for £(/) to be a Lie algebra of maximal dimension 3 we must have 
fx(x) — f x {[x, y\) — f v {y) — fy([x, y]) — 0, and we have seen before that f y (x) must be equal to 
f x {y). This means that (F, r) is not a free /-set, since ip(v) is not a Lie algebra if V( Xl!/ ) ^ V( y , x )- 
However, if we take F = {(x,y)}, r x (y) = r y (x) = i x (y), and r x (x) = r x ([x,y]) = r y (y) = 
r y (x,y) — 0, then (F, r) is a sufficient /-set of size 1, and by checking the Jacobi identity for 
£(/) we find that it is a free sufficient /-set. 

Lemma 10. IfF is a free sufficient f-set then X = K F . 

Proof Recall from above the map ip : K F — »■ (J 7 *) 11 . By the assumption that F is free, we have 
ip(K F ) C X. The assumption that F is sufficient guarantees that X is not a proper subset of 
ip(K F ), so that ip is a bijection between K F and X. □ 

We remark that a priori it is possible for X to be an affine space while no free sufficient /-set 
exists. However, we have not yet encountered such a case (see also the discussion in Section 7). 

5 The algorithms used 

In this section we describe the algorithms used to obtain the results in Section 6. For the com- 
putational results in this paper we have developed an algorithm that, given a graph T and a field 
K, computes the variety X and the multiplication tabic of £(/) for any / e X. These algo- 
rithms were implemented in the Magma computer algebra system [BC10]. The algorithm may 
be broken up into four distinct parts, of which the third (MultiplicationTable) and the fourth 
(MinimizeFSet) are the most time-consuming. 

5.1 Computing a basis for £(0) 

In this section we describe how we compute, strictly speaking, a set of monomial elements of 
J-{K, T) that projects to a basis of £(/) for any / e X. However, because in particular such 
a basis projects to a basis £(0) that consists of monomials, we will often call such a basis "a 
monomial basis of £(0)" in the remainder. 

A sketch of the algorithm ComputeBasis is given as Algorithm 1. Wc initialize B with the 
set n of generators of T(K, Y), so that B contains a basis for all monomials of length 1. Wc then 
iteratively extend a basis for monomials up to length I to a basis for monomials up to length I + 1 
by forming all products [x,b], where x € n and b a monomial of length I, and testing whether this 
product can be written in previously found basis elements. Once we have arrived at a length for 
which no new basis elements are found, we are finished. 

Note that while we return elements of J 7 , all nontrivial computations take place inside the free 
associative algebra U on the sandwich generators of a Lie algebra L (which will be identical to 
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ComputeBasis 

in: A field K and a graph T, with II = V(r). 
out: A set of monomial elements of T . 
begin 

1 let U be the free associative algebra with generators II, 

2 let fx : T — > U be the linear map defined by x i-> x for x € II and 

[6, c] i-» fi{b)n{c) - n{c)n(b) for 6, c e £(0), 

3 let B = B' = n C T and By = B' v = {(i(x) \ x G B} C [/, 

4 let C/D = (x 2 | xe n) U (/i([x,y]) | x,y e n,x^y), 

5 while ^ do 

//Update the ideal representing extremality of elements ofU 

6 let I U = I U U (n([x, [x,b]]) | x e n,6 e B'), 

//Find new elements that are not linear combinations of existing basis elements 

7 let B' Q = B' C T, B' = C T, B' v = C U, 

8 for each x e II, b e B' do 

9 let p = 6]), 

10 if p £ KBu + I v then 

n let B = BU [x,b], B' = B' U [x,b], 

12 let B U = B U U {p}, By = B[j U {p} 

13 end if, 

14 end for, 

15 end while, 

16 return B. 
end 

Algorithm 1: Computing a basis 
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InitialFSet 

in: A field K and a graph T, with II = V"(r), and a monomial basis B of C(0). 

out: A sufficient /-set (F,r). 

begin 

1 let F = 0, 

2 for c e -B, y e II do 

3 if f y (c) can be expressed in {f x (b) \ (x,b) < (y,c)} using Lemma 9(i) - (viii) then 

4 set r y (c) e Rp to be the corresponding expression, where f x (b) is replaced by r x (b) 

5 else 

6 set F = FU{(y,c)}, 

7 set r y (c) = f y (c) 

8 end if 

9 end for 

10 return (F, r). 
end 

Algorithm 2: Finding an initial /-set 

£(0)). In order to connect the two algebras we introduce the map \i in line 2, that maps elements 
of T to elements of U. In the course of the algorithm, we are also constructing an ideal Ijj of U 
such that the adjoint representation U — >• End(X ) factors through U —tU/Iu- Here ly describes 
the fact that the generators are sandwich elements; Ijj should be viewed as the analog in U of the 
ideal 1(0) of T (cf. Section 2). This construction is what allows us to include x 2 in Ijj in line 4 
of the algorithm, since ad^j = if a; is a sandwich. 

Note finally that in order to be able to compute in U/Iu (as required in line 9 of Algorithm 1) 
we have to repeatedly compute a Grobner basis of Ijj but, fortunately, since the elements of Ijj 
are all homogeneous a truncated Grobner basis is sufficient (and sufficiently efficient in the cases 
we are interested in). 

Lemma 11. The algorithm ComputeBasis returns a set B of monomials in T , and B projects 
to a basis of L(f) for any f G X . Moreover, all elements of B are either of the form x (x € nj 
or of the form [x, b] (x £ II, b e B ). 

Proof The fact that B consists of monomials, and that these monomials are of the stated form, is 
immediate from the algorithm (note the construction of new elements of B in line 11). Moreover, 
because our / £ X arc such that dim(£(/)) = dim(£(0)), a homogeneous basis for £(0) is a 
homogeneous basis for £(/) (cf. [itpPR09, Lemma 3.1].). □ 

5.2 The initial sufficient /-set 

The computation of the multiplication table, described in Section 5.3, will take place over the ring 
i?F for some sufficient /-set (F, r). In this section we describe InitialFSet, the procedure for 
initialisation of Rp, as presented in Algorithm 2. Note that the /-set returned by this algorithm 
is a sufficient /-set, but in general not a free /-set. 

We fix in advance an (arbitrary) total ordering on the elements of II, and a total ordering on 
B respecting the natural order by monomial length, i.e., b < c whenever 1(b) < 1(c). These extend 
to a total ordering of II x B by (a;, b) < (y, c) whenever b < c, or b — c and x < y. The main for 
loop in Algorithm 2 traverses II x B in ascending order; this is the ordering meant by "<" in line 
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MonomialToBasis 

in: A field K and a graph T, with U = V(T), a monomial basis B of £(0), 
a sufficient /-set (F,r), a partially defined multiplication table (a F bc ), 
and a monomial m = [x\, . . . ,x{\. 

out: w G (Rf) B such that m = X&es v ^>-> 01 f an - 

begin 

1 if m can be written as a linear combination w G (Rf) B of basis elements, 
by iteratively computing [x^, . . . , X;] for i = Z — 1, . . . , 1 using (a F fcc ) then 

2 return v. 

3 else if x;_i 7^ x; in T then 

4 return G (i? F ) S - 

5 else if = x i+1 for some i < I then 

/* Use [x,[x,y]]=f x (y)x */ 

6 let x = Xi, y = [x i+2 , . . ■ ,xi], 

7 let t = r x {y) G R F , 

8 let v such that ^ ^ = - ■ ■ (using (a F bc )), 

9 return t ■ v. 

10 else if Xi = x i+2 for some i < I — 2 then 

/* f/ S e2[x,[y,[x,z]]] = f x ([y, z])x - f x (z)[x,y] - f x (y)[x, z] */ 

11 let x = Xi, y = x i+ i, z = [x i+3 , . . .,x{\, 

12 let t xyz = r x ([x i+ i,z]) G R F , 

13 let t xy = r x (y) e -Rf, 

14 let = r a (z) G i? F , 

15 let w 2 such that v b b = [ x 1j ■ ■ ■ > x i\ ( us i n g ( a x6c))' 

16 let v xy such that X>j^6 = fci' ^2, ■ ■ • > N-i, [ x i,y]]]] (using (a F 6c )), 

17 let u* 2 such that Y, v b Zb = l x ii I x 2, ■ • ■ , [zj-i, z]]]] (using (a F bc )), 

18 return • - \t xz ■ v xy - \t xy ■ v xz . 

19 else 

20 return fail 

21 end if. 
end 

Algorithm 3: Expressing a monomial as a linear combination of basis elements 

3 as well. InitialFSet uses the relations in Lemma 9 as reduction steps: A = B is interpreted 
as " A may be reduced to B" . 

Lemma 12. The algorithm InitialFSet returns a sufficient f-set. 

Proof Observe that the relations in Lemma 9 all must hold if C(f) is to be a Lie algebra of 
dimension £(0); this immediately implies sufficiency. Finally observe that all reductions thus 
obtained are polynomial, so that indeed r y (c) defined in line 4 of Algorithm 2 is an element of Rf- 
□ 

5.3 Computing the multiplication table 

The third part of the algorithm is MultiplicationTable, where part of a multiplication table 
for £(/) is computed. Let (F, r) be a sufficient /-set (such as the one produced by InitialFSet), 
and let II be the set of generating extremal elements and B be a monomial basis of £(0) as before. 
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MonomialRelations 

in: A field K and a graph T, with II = V(T), a monomial basis B of £(0), 
a sufficient /-set (F,r), a partially defined multiplication table (a^ bc ), 
a set M. of monomials of length / currently under consideration, 
an integer k > 1, and a monomial to = [xi, . . . , x{\ € M.. 

out: A set A of elements a e (Rp) BxM such that to = J2beB ab ^ + Smga-i a " u - 

begin 

1 let A = 0, 

/* By the Jacobi identity */ 

2 write [xk, ■ ■ ■ , x{\ as a sum s of 2'~ fe+4 monomials using Equation (15), 

3 find a e (i? F ) Bx - A/1 such that [xi, [x 2 , . . . , [x fe _i, s]]] = Efces a b fo + J2ueM a « w ' 
using MonomialToBasis to obtain the entries indexed by B, 

4 let ,4 = A U a, 

/* Using Lemma 13 */ 

5 for such that Xi = Xj, 1 <i <k, and i + 3 < j < I do 

6 write to as a sum s of i — j + 1 monomials using Lemma 13, 

7 find a e (i? F ) Sx ' M such that s = X^ G s a b fe + S«eA-i a « u ' using 
MonomialToBasis for the entries indexed by B, 

8 let A = A U a 

9 end for, 

10 return A. 
end 

Algorithm 4: Finding monomial relations 
We seek to compute a^ bc <G Rf such that, for x, b e _B, for specialisation to any / e X, we have 

We will call such a set (a^ bc ) a multiplication table of £(/). 

Note that this is a multiplication table over i? F ; for / £ X it specializes to a multiplication 
table for the Lie algebra £(/) over K upon evaluation f y (c) i-> /^(c) for (y, c) € F. First, observe 
that once we have computed o^ bc for all x € II and all 6 G B every monomial [x\, . . . , x{\ can easily 
be written as a linear combination of basis elements by using (a^ bc ) to multiply "from the inside 
out" . Secondly, observe that then similarly for any (6, c) € B x B we can use the Jacobi identity 
to write [6, c] as a linear combination of monomials, and then apply the same trick (it would of 
course in general yield an exponential number of terms). Therefore we restrict to computing the 
(4) for x G n, 6 e B. 

We present sketches of the algorithms used as Algorithms 3, 4, and 5, where 5 is the main 
function. We clarify some of the notation: Firstly, we consider elements of (Rf) B as vectors 
indexed by elements of B, as for example in lines 8, 15-17 of Algorithm 3. These vectors may 
of course be multiplied by elements of i? F , as in lines 9, 12-14 of that algorithm. The same 
notation is used for elements of (Rf) M and (R F ) BxM - Secondly, we write e m for the vector with 
everywhere except for the position indexed by to (cf. line 10 of Algorithm 5), and a\s for the 
restriction of the vector a to a subset S of its index set. 

To clarify the monomial rewriting in line 6 of Algorithm 4, we state the following lemma. 



10 



ComputeMultiplicationTable 

in: A field K and a graph T, with IT = V(T), a monomial basis B of £(0), 

and a sufficient /-set (F, r). 
out: A partial multiplication table (a^ c ), where (x, 6, c) e II x B x £>. 
begin 

1 let n — maxbgs \b\, 

2 for i = 1, . . . , n do 



/* Consider monomials of length I */ 

3 let M = {[x,b] \xeT\ 1 beB such that |b| = I}, 

4 let M L be the empty x |A4| matrix, 

5 let M R be the empty x \B\ matrix, 

6 let k = 1, 

7 while rk(M L ) < \M\ do 

8 for m E Ai do 

/* -Fmd relations for m */ 

9 if v = MoNOMlALToBASls(m) ^ fail then 

10 append e m to M L , append v to M R 

11 else 

12 let A = MonomialRelations(A1, k, m), 

13 for each a e A append e m — a|x to M L , append a|s to M 

14 end if, 

15 end for, 

16 let k = k + 1, 

17 end while, 

/* Read the result */ 

18 for m = [x,b] £ M do 

19 let v = wM R , where w is such that wM L = e m , 

20 let a* bc = v c for c E B 

21 end for 



22 end for, 

23 return (o£ 6c ). 
end 

Algorithm 5: Computing a multiplication table 
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Lemma 13. Let m = [x\, . . . ,x{\ be a monomial of length I, such that Xj = xj for some i,j such 
that j — i > 3. Then m may be rewritten as a sum of 3 monomials of length less than I and j — i — 2 
monomials of length equal to I. 

For brevity and ease of reading we will prove the following easier lemma, of which the above 
is a straightforward generalisation. 

Lemma 14. Let m = [xi, . . . , x{\ be a monomial of length I, such that x\ = Xj for some j such 
that j > 4. Then m may be rewritten as a sum of 3 monomials of length less than I and j — 3 
monomials of length equal to I. 

Proof Observe that we have the following equality of a monomial of length I > 3 and a sum of 
different monomials, obtained by repeatedly applying the Jacobi identity: 

[Xi, [x 2 , [x 3 , . . . ,X;]]] = [x 2 ,[xi,[x 3 ,...,Xj]]] - [[x 3 , . . . ,X;], [xi,X 2 ]] 

= [x 2 ,[xi,[x 3 ,...,X*]]] - [x 3 ,[[x 4 ,X 5 ,...,X;],[xi,X 2 ]]] 

+ [X4, [[x 3 ,x 5 , . . . ,Xl], [xi,x 2 ]]] 
= . . . (2^ 3 terms in total). (15) 

Moreover, it follows immediately from the Jacobi identity that 

[Xi, [x 2 , [x 3 , . . . ,Xi]]] = -[[x 2 ,Xi], [x 3 , . . . ,X;]] + [x 2 , [Xi, [x 3 , . . . , X;]]] (16) 

Observe that there are i — 2 generators between x\ and Xi, so that by applying Equation (16) 
z — 3 times (but starting at x 2 rather than xi), we find 

[%i, [X2, [x3, ■ ■ -,xi}]] = [xi, [[xi-i, . . . ,x 2 ], [xi, [x i+ i, . . ■ ,x/]]]] + [i - 3 monomials of length I], 

(17) 

By Equation (6) the first monomial reduces to monomials of smaller length (recall Xi = Xj): 

[Xi, [[Xi_i, . . . , X 2 ], [Xi, [X i+1 , . . . , X/]]]] = i/ Xl ([Xj_i, . . . , X 2 ], [Xj+i, . . . , Xl])X! 

- \ fxA\Xi+l, ■ ■ -,xi])[xi, [Xi-1,.. • ,x 2 ]] 

- ^fxAiXi-1, ■ ■ ■ ,X2])[xi,[x i+1 , . . . ,Xi\], (18) 

so that, combining Equations (17) and (18), we find 

[xi, [x 2 , [x 3 , . . . , x;]]] = [3 monomials of length < I] + [i — 3 monomials of length /], (19) 

proving the lemma. □ 

Lemma 20. //MultiplicationTable returns a multiplication table (a^ bc ) then, for any f e X, 
these {a^ bc ) specialize to a multiplication table for the Lie algebra £(/) upon evaluation f y (c) 
f y (c) for (y,c)eF. 

Proof Recall that, as described at the start of Section 5.3, the structure constants (a^ bc ) for 
x, b, c G B follow immediately from those where x € n and b,c e B. We will consequently assume 
that (a^ bc ) is known for x,b,c G B. 

Observe that in Algorithm 5 the products are computed by ascending length, and for each 
length I two matrices M L and M R are computed. Throughout the algorithm the rows of M L 
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MinimizeFSet 

in: A field K and a graph T, with II = V(T), a monomial basis B of £(0), 

a sufficient /-set (F,r), and a multiplication table (a£ 6c ). 
out: An updated /-set (F,r) of smaller size, and updated multiplication table, 
begin 

/* Collect relations between elements o/F */ 

1 let ft = {[y, [y, c]] - f y (c)y \ (y, c) e P}, 

2 let ft = ft U {[a, [6, c]] + [b, [c, a]] + [a, [b, c]] \a,b,ce B}, 
/* Reduce size of F using ft */ 

3 do 

4 for each non-zero coefficient t E Rf of each element u <G ft do 

5 find a e if*, and (y,c) such that f y (c) occurs only in the linear term af y (c) of t 

6 or continue if no such element exists, 

7 setr v (c) = -±(t-af y (c)), 

8 replace i y {c) by r y (c) in ft, 

9 replace f y (c) by r y (c) in (a£ 6c ), 

10 set F = F\{(y,c)} 

11 end for 

12 until F is unchanged, 

13 return (F,r), (a^). 
end 

Algorithm 6: Finding a minimal sufficient /-set 

and M K contain linear relations between elements of M and elements of B, i.e., the i-th row 
of M L must be equal to the i-th row of M R . To see that these relations are true, observe that 
those gained from MonomialToBasis (Algorithm 3) are true by induction on the length of m 
(if returned in line 2), by definition of extremality (if returned in line 9), or by the first Premet 
identity Equation 6 (if returned in line 18). Similarly, those gained from MonomialRelations 
(Algorithm 4) are true by the Jacobi identity (those produced in line 3) or by Lemma 13 (those 
produced in line 8). 

Also, because we traverse the products of basis elements by ascending length, the elements 
of (Rf) B required in lines 8 and 15-17 of Algorithm 3 are guaranteed to be easily deduced from 
the entries of (a^ 6c ) computed earlier. Finally correctness of the a^ bc computed by solving the 
system of linear equations represented by M L and M R in line 18-21 of Algorithm 5 follows from 
bilinearity of the Lie algebra multiplication. □ 

We remark that we cannot give an estimate of the number of relations we have to produce 
(i.e., the number of times we call MonomialRelations and the value of k) before M L is of 
full rank. In practice, however, we find that k docs generally not exceed 2. In particular, this 
implies that termination of ComputeMultiplicationTable is not guaranteed - but by Lemma 
20 correctness upon termination is. 

5.4 Finding a minimal sufficient /-set 

After the multiplication table has been computed, we try to minimize the size of the /-set. This 
procedure is called MinimizeFSet and presented here as Algorithm 6. The procedure is straight- 
forward: we find a set of relations ft between the f y (c) by requiring that the Jacobi identity holds 
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X: 

dim(L): 
Lj Rad(L): 



K 1 

3 



3-dim: 
Os 



runtime: 




X: 

dim(L): 
Lj Rad(L): 



K 2 
G 



(c) 




X: 

dim(L): 
L/Ra,d(L) 



8 



3-dim: 
Os 



A: 



8-dim: A 2 
Os 



runtime: 



runtime: 



Table 1 



Computational results (2 or 3 generators) 



and that for all y € II, c € B we have [y, [y, c]] = f y {c)y, and we subsequently use those relations 
to find nontrivial r y (c) and remove (y,c) from F. 

Once the algorithm finishes, we apparently could not find any new linear relations. This does 
not necessarily mean that the /-set (F,r) returned is free. However, if it happens to be the case 
that 7Z — {0} upon exiting, then (F,r) is indeed a free sufficient /-set since then the Jacobi 
identity apparently holds. For all the cases we have tested, the /-set returned was in fact free. 

Observe that in this algorithm only division by elements of K* is taking place (cf. lines 5- 
6 of Algorithm 6) and not by other elements of Rf, so that all these calculations take place in 
multivariate polynomial rings Rf over K. While it is a priori conceivable that division by elements 
of Rp\K* would be necessary to reduce the size of the /-set (and thus that our algorithm would 
return a non-free /-set), we have not observed such a case in practice. 

We remark that a rather more sophisticated alternative to Algorithm 6 has been implemented. 
For example, we execute a similar algorithm in the course of the computation of the multiplication 
table (rather than only after the fact), wc alternate between collecting a relatively small number 
relations and reducing the size of F (rather than do them in succession), and in line 5 of Algorithm 
6 we give preference to y, c such that the pair (y, c) is large in the ordering on II x B introduced 
in Section 5.2. 

Lemma 21. IflZ = {0} upon exiting of the algorithm MinimizeFSet, then X is isomorphic to 
the affine space K' F \ . 

Proof This follows immediately from Lemma 10. □ 

6 Computational results 

We present the results of applying the algorithms described in Section 5 in Tables 1-3. For 
each connected undirected finite graph T with n vertices, 2 < n < 5, without loops or multiple 
bonds, we state the dimension of the variety X, the dimension of £(0), the generic dimension and 
isomorphism type of Lj Rad(L), and the CPU time taken (on a 2 GHz AMD Opteron) to execute 
the steps described in Section 5. 

Recall that we claim the results presented in Tables 1 - 3 to be valid for any field K of 
characteristic distinct from 2. To establish this, we have performed the following calculations for 
every graph T occurring in the results. First we calculate a basis of the equivalent Lie ring over 
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A: 

dim(L): 
Lj Rad(L): 
runtime: 



A 3 
12 

3-dim: Ai 
Os 



(d) 




A: 

dim(_L): 

L/Rad(L): 

runtime: 



A 5 
15 

15-dim: A 
Os 



•- 



X: 

dim(i): 

L/Rad(i): 

runtime: 



A 3 
10 

10-dim: B 2 
Os 



(e) 




A: 

dim(L): 
L/ Rad(L): 
runtime: 



A 8 
21 

21-dim: B 
Os 



A: 

dim(L): 
L/ Rad(L): 
runtime: 



A 5 
15 

15-dim: A 3 
Os 



(f) 




A: 

dim(L): 

L/Rad(L): 

runtime: 



A 12 
28 

28-dim: D 
Os 



Table 2: Computational results (4 generators) 



Z using a procedure similar to the one described in [CdG09, Section 3]. This gives us a set V\ of 
primes for which the monomial basis computed in Algorithm 1 would differ from the characteristic 
case. Then we execute the algorithms described in Section 5 for K = Q, and in Algorithm 6 we 
store in the set Vi the prime components of the denominator of the a we divide by. Subsequently, 
we perform the computation with A = GF(p) for p e (Pi U 7 , 2)\{2}; fortunately this last task 
is easily done in parallel. Only in one case did we find a different result: for T = K 5 the case 
where char(A) = 3 is of dimension 538 (rather than 537). The result that A = {0} remains valid, 
however. 

We emphasize that in each of the cases considered we found a free /-set, and consequently 
prove that A is an affine space (cf. Lemma 10). Moreover, we keep track of the reductions and 
relations used to turn the trivial sufficient /-set into a free sufficient /-set, so that we obtain a 
certificate for the correctness of these results. Unfortunately, however, such certificates are far too 
long to reproduce here. 

Before turning to some of the results, let us elaborate on a particular piece of data, namely the 
isomorphism type of Lj Rad(L). For a number of special T (namely the Dynkin diagrams of affine 
type) it has been shown that the Lie algebra £(/) is isomorphic to a fixed Lie algebra L for / in an 
open dense subset of A [Ditp08, Theorem 22]. A similar result is proved in [itpPR09, Section 7] 
for the four infinite families there considered. We investigate this property in our cases as follows. 
Firstly, we carry out all the computations described in the previous section over Q. Recall that in 
all cases we found a free /-set (F,r); we let I = |F|, so that A = Q 1 . Now take some finite field 
F g , whose characteristic p is distinct from 2 and such that no multiples of p occur in denominators 
in the multiplication table (a^ bc ) (in our examples the multiplication table always turns out to be 
integral, so any finite field of characteristic distinct from 2 works). Then take v uniformly random 
in ¥ l q and let L — £(ip(v)), where we interpret the multiplication table {a^ bc ) over Rp (g> ¥ q and 
is evaluation as described in Section 4. This gives us a Lie algebra L defined over F 9 , for which 
it is straightforward for Magma to compute the isomorphism type of L/Rad(L). We repeated 
this process a number of times over various finite fields (e.g., F 17 , F 10 i) and stored the results - in 
each case finding the same dimension and isomorphism type of L/Rad(L) in the vast majority of 
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X: 

dim(L): 
Lj Rad(L): 
runtime: 



X: 

dim(L): 
Lj Rad(L): 
runtime: 



K 5 
28 

28-dim: D 4 
Os 



K 4 
20 

10-dim: B 2 
0s 



(g) 



(h) 



6 



X: 

dim(L): 

L/Rad(L): 

runtime: 



X: 

dim(L): 

L/Rad(L): 

runtime: 



X 6 
30 

15-dim: A 3 
0s 



K 6 
24 

24-dim: A 4 
0s 



X: 

dim(L): 
L/ Rad(L): 
runtime: 



15 

10-dim: 
0s 



B 2 




X: 

dim(L): 

L/Rad(X): 

runtime: 



52 

52-dim: F 4 
0s 



X: 

dim(L): 
Lj Rad(L): 
runtime: 



K 7 
36 

36-dim: 
0s 



B, 



0) 




X: 

dim(L): 
Lj Rad(L): 
runtime: 



K 9 
45 

45-dim: 
0s 



D, 



X: 

dim(i): 
Lj Rad(L): 
runtime: 



K 6 
30 

15-dim: A3 
0s 



(k) 




X: 

dim(L): 
Lj Rad(L): 
runtime: 



K 9 
45 

45-dim: D 5 
0s 



X: 

dim(L): 
Lj Rad(L): 
runtime: 



24 

24-dim: 
0s 



A 4 



(1) 




X: 

dim(L): 

L/Rad(X): 

runtime: 



52 

52-dim: F 4 
0s 



Table 3: Computational results (5 generators) (1/2) 
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X: 

dim(L): 
Lj Rad(L): 
runtime: 



K 9 
45 

45-dim: D 5 
Os 



(r) 




X: 

dim(L): 

L/Rad(L): 

runtime: 



K xi 
134 

28-dim: D 4 
10s 



X: 

dim(L): 
Lj Rad(L): 
runtime: 



K 13 

86 

28-dim: D 4 
3s 



(s) 




X: 

dim(L): 

L/Rad(X): 

runtime: 



K 21 

133 

133-dim: E 7 
14s 



X: 

dim(i): 
Lj Rad(L): 
runtime: 



78 

78-dim: E 6 

2s 



(t) 




X: 

dim(L): 
Lj Rad(L): 
runtime: 



K ix 
249 

78-dim: E 6 
2510s 



X: 

dim(L): 
Lj Rad(L): 
runtime: 



X: 

dim(L): 
Lj Rad(L): 
runtime: 



K 14 
78 

78-dim: E 6 
2s 



78 

78-dim: E 6 
2s 



(u) 




X: {0} 

dim(i) • I 538 if char(X) 
[ 537 otherwise 

L/Rad(i): trivial 

runtime: 38260s 



Table 3: Computational results (5 generators) (2/2) 
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cases. We remark that we have not distinguished between split and twisted forms of Lie algebras 
of the same Cartan type. 

The following results are particularly worth pointing out (we refer to the case presented in 
Table i(x) simply as case [ix]). 

(i) [lc], [2c], [3a], and [3h] agree with cases proved in [Ditp08]; 

(ii) [2b], [2d], [3f], and [3j] agree with cases proved in [itpPR09]; 

(iii) [3n]-[3q] are the four cases with 5 vertices and 3 commuting edges. Where one might expect 
a structural similarity, [3n] is the odd one out both in terms of dimension (86 vs 78) and 
type (a large radical vs simple of type E 6 ); 

(iv) Similarly, [3r] and [3s] are the two cases with 5 vertices and 2 commuting edges. Again, 
despite the difference in dim(£(0)) being only 1, their isomorphism types are very different; 

(v) In our computations, we have not encountered a Lie algebra of type Eg. Even though case 
[3t] is of dimension 249, the biggest simple component we found there was of type E 6 . It 
has been proved that Eg can be generated by 5 extremal elements [CSUW01, Theorem 8.2], 
but apparently not in the manner considered in this paper: Eg is not a subquotient of 
any maximal-dimensional Lie algebra generated by 5 extremal elements plus commutation 
relations, for then it would have to occur either in [3t], where it does not, or in [3u], where it 
cannot for such a Lie algebra is always nilpotent. On the other hand, Eg can be generated 
in this manner by 9 extremal elements [Ditp08, Theorem 22]. 

(vi) [3u] is the case of largest dimension among the 5 generators; it is also the only case where 
X is a point and consequently only one Lie algebra occurs, and that it is a nilpotent Lie 
algebra (cf. [CSUW01, Lemma 4.2]). 

7 Conclusion 

Proof of Theorems 3 and 5. These theorems follow immediately from the correctness of our 
algorithm, as proved in Lemmas 11, 12, 20, and 21, and the computational results presented in 
Section 6. □ 

In Theorem 3 and Corollary 4 we prove that if T is the complete graph on 5 vertices, then 
X = {0} and the only Lie algebra that occurs is nilpotent. This of course leads to the question 
whether this pattern continues: "For F the complete graph on n vertices (n > 5), is X a point?". 
Unfortunately, this question is not feasible using our software (experiments show that for n = 6 
already dim(£(0)) > 140000). Another approach would be to consider relations between the 
variety for T and the varieties for subgraphs F' C T. This problem is the subject of ongoing 
research. 

Secondly, as mentioned earlier, two questions were posed in [Ditp08, Section 5.2]: "Is X always 
an afhne space?" , and "Is there always a generic Lie algebra?" . Although the authors of that paper 
expect the answers to both questions to be negative, we have found no counterexamples. Moreover, 
since these questions can be answered affirmatively for both a number of infinite series of sparse 
graphs (as studied in [Ditp08, itpPR09]) and a number of dense graphs (as studied here), we would 
not dispose of these conjectures too easily. 
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